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Abstract 

The possibility that a static magnetic field may decay through 
production of electron positron pairs is studied. The conclusion 
is that this decay cannot happen through production of single 
pairs, as in the electric case, but only through the production of 
a many-body state, since the mutual magnetic interactions of the 
created pairs play a relevant role. 

The investigation is made in view of the proposed existence 
of huge magnetic field strengths around some kind of neutron 
stars. 



* E-mail: giorgio@ts.infn.it 



1 



1. Introduction 

The perturbative calculations in QED have been pushed to an exceptional de- 
gree of refinement as far as their comparison with the high precision experiments. 
The study of the situations which require a genuine non-perturbative treatment is 
less systematic, also in this field, however, there are problems that can be consid- 
ered classical and have gained a high degree of clarification. One of these prob- 
lems concerns the effects of intense external field: the decay of the vacuum subject 
to a strong electric field through emission of e~ e + pairs is well settled through 
the analysis given by Schwinger, starting from the Euler-Heisenberg effective La- 
grangian[l,2]. If one considers strong magnetic fields there is a strict analogy if 
we admit the existence of monopoles, in that case the problem is only quantitative 
since the monopoles are expected to be very heavy, but for ordinary particles the 
parallel with the electric case is less straightforward. 

The difference is that we associate a potential energy to the charge in an electric 
field, but we cannot do the same for an electric charge in a static magnetic field. 
The possible instability of the vacuum in the magnetic field, on the other side, seems 
interesting to be analyzed in view of the guess that extremely high magnetic field, of 
macroscopic extension, are realized in nature, around some neutron stars [3]. If the 
magnetic field is not completely static there is certainly a pair production, it can be 
seen e.g. treating the problem with the formalism of the adiabatic approximations^], 
but in this case the rate of production depends on the square of the time derivative 
of the field. Here a qualitative and partially quantitative analysis of the possibility 
that the some kinds of vacuum instability leads to pair production even in static 
conditions is presented. 

The main idea is stated in this way: any number of pairs in a given magnetic 
field cannot give rise to instabilities unless we take to some extent into account 
also their mutual interaction: in fact in this last case every pair shields partially 
the magnetic field in which the other particles lie, the overall effect could be a 
decreasing of the field intensity. This effect, when the original field is very strong, 
may decrease the density of magnetic energy and compensate the cost in energy 
for the creation of the pairs of charged particles. If this is true the effect is largely 
collective and cannot be studied particle by particle. A simple energy balance of 
this collective production is presented in section 2 and the reason for believing in the 
instability are given, since at this point we are in presence of an e~ e + multiparticle 
state the superimposed Coulomb effects are studied in section 3, some conclusions 
are presented in the last section. 

2. Main features of the model 

The physical model is described in these terms: there is a classical magnetic 
field B constant in time and in it there is a second quantized electron field The 
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total energy of the system is 



£ T = ±jB 2 d 3 r + H with H = J{^[a- (p + eA) + (3m]^}d 3 r . (1) 

The usual variational principle, for the stationary case, S£t = yields the 
equation for the magnetic field. This equation is then brought to numerical form 
by putting it among eigenstates of Ti. Two known relations are used : 

< a\d s H\a >= d s < a\H\a >= d s E a (s) (2a) 

when H depends on the parameter s, \a >is one of its eigenstates and E a is the 

— * — * 

corresponding eigenvalue. Since E will be expressed in terms of B = curlA, one 
must remember that: 

5 d 

. 7zr = e ik id k . (26) 

5Ai(r) oBiif) 

The standard relation curlB = J, yields then 

S S 
eikid k Bi = -e ik id k —E i.e. Bi = ———E + diU . (3) 
OBi OBi 

The functional derivative is zero in the limit e — > 0, so d[U represents the field in 



? (o) 
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absence of vacuum polarization, it will be denoted by B\ 

Now we specialize to a definite family of field configurations: let the magnetic 
field be uniform B = Bn and the volume of quantization of the system be a prism 
of height L, parallel to n, and square section of side R. The multiparticle state 
is obtained in the most trivial way i.e. by filling the one-particle states until some 
Fermi level and leaving completely empty the higher one-particle states. It turns out 
natural to consider two independent parameters in order to fix the highest populated 
level one for the longitudinal and one for the transverse degrees of freedom; so it 
results (see e.g. [6]): 

y = — dk-eB^-.yy 

^ 2tt J_ k 2tt ^ ^ 

levels o s=±l 

Looking for a uniform solution B^ must be a constant vector and —5E/SBi 
is also a constant playing the role of a total magnetization Mi. 



* A clear although artificial way of dealing with this kind of boundary condition 
may be found in the formalism of space-dependent coupling constants as proposed 



by Bogoliubov and Shirkov[5]. 
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The energy coming from H has a constant density e = E/V , V = R 2 L. For 
notational simplicity all the vectors are projected onto the direction n and only 
these components are now used. The total energy density is 

e T = \B 2 + e = \{B^ + M) 2 + e , (4) 

and if it can become smaller that the energy density \B^ 2 , we may argue that the 
vacuum polarization leads to an instability of the magnetic vacuum. 

The single-particle energy levels are the standard relativistic Landau levels 

w n , a {k) = \/m 2 + k 2 + eB(2n + 1 + s) 

so the eigenvalue of H is 

i r K N 

E(K,N) = 2- —LR 2 eB dk \2 ^ V™ 2 + k 2 + 2eBn + ^m 2 + k 2 ^ . (5) 

The factor 2 in front of the whole expression arises when the contribution of the 
positron states is added to the contribution of the electron states. In order to 
perform the discrete sum the approximation used is 

N o 

Ev^TG^— [F + (iV+I)G] 3 / 2 , 

n=l 

which is valid up to terms constant in N. 

In the actual situation F = m 2 + k 2 , G = 2eB. The integration over k may 
be carried out completely in a straightforward way; the result is complicated and 
not very transparent. Since we are looking for a possible instability we are free to 
choose the trial state and so the range of the parameters iV and K, a choice that 
seems promising is 

NeB » K 2 » eB , K » m , 
because it allows an expansion in decreasing powers of ^/N^B. 

e = [2K(2NeB) 3/2 + K s V2NeB] + ■■■ = aB 3 / 2 + cVB + d + ■ ■ ■ . (6) 

The coefficients a , c are defined in the expression, d denotes the addendum which 
does not contains the factor NeB. 

This expression is now used to calculate M: 

M = — ^aVB — \cj\fB~ + • • • • (7) 
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Then it is inserted into the expression for the total energy difference 



A = \B 2 + e - \[B^} 2 = e + MB - \M 2 , 



with the result: 



A 



-\aB^ 2 
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a 2 B+\cVB-{l 



ac — d) . 



(8) 



The coefficient of the leading power in B is negative and we may therefore conclude 
that there must exist configurations in which the creation of a collective state of 
electrons and positrons has the effect of lowering the total energy notwithstanding 
the cost in energy of the mass and kinetic terms for the charged particles. 

3. Effect of the Coulomb interaction 

Since we are now considering the creation of a plasma of e + ,e~ we are also 
led to consider the possible effect of the Coulomb interaction. It may be estimated 
perturbatively by inserting the corresponding two-body operator: 



between the original states. This interaction is expected to give a negative con- 
tribution to the energy and since it does not contain A it does not modify the 
magnetization. The procedure to deal with this problem can be found in standard 
textbooks [7], the actual calculations become, however, very laborious if we want 
to use the correct wave functions of the electron, which are essentially harmonic- 
oscillator functions. A simplified investigation is here presented: the particle states 
are simply represented by plane waves, the sums are cut at the same levels as in the 
previous case, so a maximum longitudinal momentum K is used and a maximum 
transverse momentum P is introduced, with the later identification P 2 = NeB. In 
this way it is possible to get a simple estimate of the Coulomb effect, in consid- 
ering the pair interaction we obtain the overall cancellation of the diagonal terms 
(e~,e~), (e + ,e + ) with (e~,e + ), which in this case is particularly evident, whereas 
the exchange terms survive and the contributions of the negative and positive charge 
add up; the result it obviously definite negative and may be expressed as: 



the last factor 2 comes here also from the sum over the two charge states. 

The integrand depends only on the difference of the momenta, so three inte- 
grations are easily performed; defining p = K 2 /P 2 the result is 






(10) 



I = 




(pi -P2) 2 



1 



d 3 Pl d s p 2 = 128K 2 P 2 



1 (l-u)(l-u)(l-«;) 
u 2 + v 2 + pw 2 



dudvdw 
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In the limit p — > Oi.e. P » K the integral develops a logarithmic singularity which 
gives the leading term: 



With the identification P 2 = NeB the leading term of the Coulomb energy is: 



Comparing it with the results of eq. (6) we see that the Coulomb energy is defini- 
tively sub-leading. This results gives us confidence that higher order perturbative 
corrections, which may also affect the magnetization, will not destroy the main 
result of the previous chapter. 

4. Conclusions 

The main conclusion we may draw from this particularly simple model is that 
the indications of an instability of the magnetic vacuum are confirmed. The source 
of this instability is to be found in the dependence of the energy density of the 
electron field on the magnetic field. The leading term grows as e Q oc B a and is 
obviously positive, the diamagnetic term has a leading addendum B(de D /dB), it 
carries a minus sign, and since 1 < a < 2 this term over-compensates the positive 
amount of energy required for the creation of pairs, the applied magnetic field 
B^ is reduced by the pair creations, it results in fact, for the actual case a = § 

B ~ B^ — |aV-B(°) + |a 2 ; the effects of the Coulomb interaction are seen only 
in the sub-leading terms. The requirement of very large magnetic field is essential 
in order that the classification in leading and sub-leading terms be meaningful, in 
fact looking blindly at eq (8) one would get the impression that the total energy 
continuously decreases with increasing N, this is clearly impossible, for a given 
B(°) at a certain value of N, the resulting field B becomes too small for the whole 
treatment to be correct. The determination of the actual values of B^ for which the 
process of spontaneous pair creation can happen is not possible within the present 
treatment; the main problem in performing a quantitative estimate of the process 
lies in the fact that it proceeds necessarily via a tunnel effect, because the creation 
of a small number of pairs is not enough to lower the magnetic energy. 

In fact the possible existence of the instability of the magnetic vacuum is strictly 
related to the mutual interaction of the created electron pairs, expressed by the fact 
that they are located in the field modified by the existence of the other pairs. The 
existence of an instability when only non interacting pairs are considered is excluded, 
in a static field, because under these conditions, the effective Euler-Heisenberg La- 
grangian, which describes the effect of virtual charged spinors in a given external 
field, never develops an imaginary part, contrary to what happens for the electro- 
static case[l,8]. The situation is also different from what expected in the spin-one 



I « 16nK 2 P 2 [\nP 2 /K 2 + const] 
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case, where for very intense magnetic fields instabilities are expected already at the 
level of single pair creation[9], but the known charged particles with spin one are 
much heavier than the electron so there should not be sizable interference between 
the two processes at the foreseen field strengths [3] of the order of 10 10 T. 

It must be, finally, noted the all the virtual effects have been ignored, they 
are likely to be important for very large field strengths, because they renormalize 
the electron charge and may also destroy the spherical symmetry of the Coulomb 
interaction[10], also the photon degrees of freedom coupled with the electron field 
may have a role, here only static interactions were analyzed, but all these and 
possibly other dynamical features are superimposed complications, they do not give, 
however, any indication of being able to destroy the main conclusion that emerges 
from the simple analysis presented in section 2. 
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